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Abstract 

The leading perturbative contributions into the two-point gauge-invariant 
correlator Tr{gF^^{x)U{x,0)U{0,y)gFp„{y)U{y,0)U{0,x)) are calculated at 
the one-loop order. It is shown, that nonlocal condensate Tr{asF^yF^y) is 
nonzero at this order. The relation with the renormalization properties of 
Wilson loops is discussed. 
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Our present theoretical understanding of strong interactions is based on 
quantum chromodynamics (QCD). It is well known, that the gauge nature of 
this theory leads to a wealth of its properties as well as to the great difficulties 
in the description of its dynamics. One source of these difficulties lies in 
the fact, that despite all observables one can calculate in the theory must 
be gauge invariant, the calculations actually involve noninvariant quantities 
like propagators at intermediate steps. The situation is studied well for the 
covariant gauges, at the same time for the noncovariant ones a lot of work 
still has to be done (see and references therein). There is also a kind of 
temptation to consider one set of gauge conditions as "more physical" than 
another one supported by real simplifications arising then one chooses the 
"most adequate gauge" for the given problem (despite all physical observables 
are of course gauge independent). 

Therefore it is very useful if possible to reformulate the theory in terms of 
gauge invariant quantities from the beginning. The latter are usually taken 
as Wilson loop functionals: 

W{C) = {Tr Pexp [ig j A^dx^)) (1) 

c 

The ordinary brackets denote usual average: 

{0[A{x)]) = J DA{x)0[A{x)]exp (-^ J d'xF^^^F^^^^ 

where field strength tensor is F^^ = d^A^ — d^A'^^ + gf°'^'^A^^Al 

Unfortunately, the loop approach has two undesirable features. First of 
all, the original equations for loop variables are written for nonrenormalized 
quantities and complete renormalization program for them seems to be very 
difficult to proceed. Second, the loop equations are complicated since they 
defined not in the familiar coordinate space but in the loop space - subtle 
mathematical object itself. 
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All that makes attempts to look for an alternative set of gauge invariant 
quantities attractive. Recently a lot of work has been done (see and 
references therein) in order to represent QCD on a way, where dynamical 
degrees of freedom are gauge invariant cumulants of the following type: 

D^ = {{f[Gix„Xo))) (2) 

G{xi,xo) = U{xQ,Xi)gF^^^X^i)U{xi,xo) 
where nonabelian phase factors are 



y 

U{x,y) = Pexp{-ig J A^{z)dz^ 



y t 



= ^-W j Mz)dz^, + {-igf J dU J dz^A^{z)A^{t) + ... (3) 

X XX 

The double brackets in denote the irreducible cumulants, for example 
{{AB)) = {AB) - {A){B), see [| and g for details. 

Applying the nonabelian Stokes theorem ^ and cluster expansion prop- 
erty 1^ one obtains: 

W{C) = {Tr Pexp {ig J A^dxf")) = {Tr Pexp {i J G{x^, Xo)da{xi))) = 

C=dS S 
oo 

= Tr Pexp {i ^ Dnd'^a) (4) 

n=l 

where Dn are defined in (0). It may be proven, that the r.h.s. of does 
not depend on the reference point Xq as it must be, if all D„ are taken into 
account. Therefore there is a natural question, which n are essential in the 
expansion (^). One easy notes, that for small contours, where the coupling 
constant is small, highest cumulants which are by definition non-reducible 
Green's functions are small too. Lattice calculations and some theoretical 
considerations demonstrate, that the same conclusion is going to hold in the 
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nonperturbative regime and only a few lowest cumulants are important. 
In particular, the lowest nontrivial bilocal cumulant measured on the lattice 
in is believed to be dominant and it is actually the case in the so called 
Gaussian model of the gluodynamics vacuum 0, §]. 

Therefore it would be interesting to understand the perturbative be- 
haviour of the lowest gauge- invariant functions (0). First of all, one could 
compare the renormalization properties of @ and (||). It will be seen below, 
that they are different, inspite of the similarities encoded in (Q). Second, such 
results may be useful for the lattice calculations, where perturbative part of 
the cumulants is measured at the small distances. It should also be noticed, 
that gluon propagator in the radial gauge is expressed as an integral of the 
quantities like (H) So analysing perturbative expansion of @ we get some 
information about the perturbative series in Fock-Schwinger (radial) gauge. 

In the present letter we are going to present leading perturbative contri- 
butions to the simplest n = 2 cumulant at the one-loop level: 

D,.p.{x, y) = Tr{gF^,{x)U{x, 0)f/(0, y)gF,^{y)U{y, 0)t/(0, x)) (5) 

The integration in any phase factor is along the straight line, connecting its 
end points. 

The detailed full calculation of this quantity in the limiting case y = 
has been done recently in ||^ by two methods - direct calculation of the 
corresponding Feynman diagrams and by methods of heavy quark effective 
theory. If the reference point [xq = in our case) is chosen at one of the 
correlator arguments (i.e. points where the operators F act are connected 
by the straight line), the quantity (|^) depends on the only vector z = {y — 
xq) — (x — xo) and two different tensor structures can be introduced: 

5f,p5ua - SupS^a = A^pa (6) 
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d 



A (1) _ ^ / f^fp s: _ , f^lf^LA _ 'I - A(2) (7) 

'-^uupcT o I 2 '^'^ 2 A*"" ^ 2 2 '^/^ / ~ '-^uvpcj \ ' J 



and 



'pupa 2 \ z"^ z"^ z"^ z"^ J f^'^P^ 

Note, that 

A^^) S S =0 

therefore only the part proportional to A*^^^ contributing to the condensate 
{agF^yF^y) . Two scalar coefficient functions in front of A*^^^ and A'^^); 

Dpupaix -y) = D{z)^^l^ + D,{z)^%^ (8) 

were actually calculated in [^j and appeared to be nonzero both contrary to 
the tree level results where only Di{z) presents. 

Generally tensor structure of (|^) is more complicated and the results of 
Q should be recovered by taking the points x, y and a;o = lying on a line, 
for example x = —y = z/2. So our results can be used as an independent 
check of §. 

Let us briefly remind the situation in QED. The abelian analog of (^) was 
investigated in [0. Calculation there simplifies greatly since phase factors 
are cancelled. Representing renormalized photon propagator in Feynman 
gauge as (we use Euclidean metric throughout the paper): 

{eA,{x)eA^{y)) = ^,.^,^ (9) 

where z = y — x and d{z'^) = e^{z^) - the renormalized effective charge and 
differentiating (R|) with respect to z one easily finds: 



{eF^,^{x)eFp„{y)) = ^ 



A (2) (^_^^^\ _ A (1) ^ 



(10) 



Here z should be taken nonzero in order to avoid contact terms and the 
derivative d' = de^{z^)/dz^. Only the term proportional to d" contributes to 



(Ff^iyFi^iy). Using the definition from |jTO| for the coordinate /3-function 



it is easy to see that 

{e'F,.{x)FM) = ^ ^ (l - ^ ) ■ (H) 



Taking e^{z^) = e^ijj.^) + be'^{iJ^)ln[l/{fizy] + O(e^) we find at the lowest 
nontrivial order: 

{e'F,^{x)F,^{y)) = -^^^ (12) 

At the tree level in the SU (N) Yang-Mills theory the only difference from 
abelian case in the common colour factor (A^^ _ ^^2: 



r){0) ^ ~ 1 9^ ]_ . A{2) QON 



Situation at the loop level is much more complicated in the nonabelian case. 
The main reason is the presence of phase factors U{x, y) in the expression (^. 
Integrations over the straight lines, connecting the points x, y and xq = lead 
to a new type of divergencies, additional to the usual ultraviolet ones. These 
divergencies are studied well in context of renormalization of the Wilson 
loops and we will explain below the role they play in renormalization of field 
correlators. 

Let us proceed with the direct calculation of the leading g'^/e terms in (^. 
The Feynman diagrams, contributing to the quantity under consideration are 
presented in Fig.l and Fig. 2. Naively only diagrams la, b, d and 2a, 6, c, /, h 
give the leading contribution. Tree level contribution D^^J^^ is defined in (pISl). 
Using dimensional regularization d = 4 — 2e and keeping only 1/e terms we 
get in the Feynman gauge (we consider gluodynamics but one can include 
dynamical quarks easily): 
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2b) + 2h) = '^N 



^d) = ^N Q)4[(^-7)ct^?(7)]-4°U 
In the last expression 7 is an angle between x and y. It may be seen that 
diagram 2a does not contain singular parts and diagram 2/ equals to zero. 
The straightforward check shows the coincidence with the results for the 
leading logarithmic terms from 0. Note, that contrary to we are working 
in the fundamental representation and our D from (^) is in fact V/2 from 

i- 

All the above diagrams proportional to A*^^^ and therefore do not con- 
tribute into the condensate 

it is fullfilled by the only 2c, the result is: 

D i^iu pa (^-^ J y) ' ^pp^ua 



-N{N' - 1) - ■ - 



A- 



x'^z^ 



+ 



\_^{yz) 



2n 



(14) 



167r4 ' ' \t) z^ 

We note the nontrivial dependence on the angles between x, y and z = y — x, 
expressed via scalar products (xz) and (yz). In particular case, if the special 
position of the points x, y and is chosen the leading logarithmic contribution 
to the condensate is absent and only finite renormalization is needed. Namely, 
it is the case if the points x, y and form the equal-side triangle of arbitrary 
size. This circumstance might be useful in the lattice calculations for fixing 
the normalization scale of the leading logarithms. 

We see that the one-loop contribution to the field correlators cannot be 
entirely explained by charge renormalization. It is clearly indicated also in 
1^ where besides terms proportianal to the first coefficient of the /5-function 
there are additional terms in Di and D. 
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To proceed, let us compare the renormalization procedure for the Wilson 
loops (see |11] for details) and correlators. The former depends on whether 
the corresponding Wilson contour C is smooth or not despite the final results 
look in some sense similar. In the former case the average is renormalizable 
in the following sense: 

W{C) = Z ■ WrUC) (15) 

where infinite Z-factor contains only linear divergencies arising from the in- 
tegrations over the contour: 

Z~exp|c J dx^ J dyuD^^,{x - y)\ ^ exp (^c^^ (16) 
\ c c / 

where L is the length of the contour a - ultraviolet cutoff and c - numerical 
constant not important for us here. At the same time logarithmic divergen- 
cies are absorbed in a usual way, so that renormalized Wilson average Wren in 
the Eq. (|T5|) is functional of the contour C and function of the renormalized 
coupling constant grenifJ') defined on the corresponding dynamical scale /i. 

It is worth noting that a single diagram in the perturbative expansion of 
the r.h.s. of (|T5|) typically gives contributions to both Z-factor and Wren{C), 
i.e. contain linearly and logarithmically divergent parts. Hence the factor- 
ization in the r.h.s. of ( P^ ) has double-faced nature - it is factorization of the 
diagrams together with the separation of contributions, coming from dan- 
gerous integration regions (i.e. linearly divergent terms) for each diagram 
itself. 

As a result contour divergencies can be separated as a common Z-factor, 
which physically renormalizes the bare mass of the test particle, moving 
along the Wilson contour. This additional infinite constant in the exponent 
reflects the impossibility to observe the Wilson loop itself. At the same time 
if one extracts the physical quantities from the Wilson loop average then only 
logarithmic renormalization cause observable effect. 
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If the contour possesses cusps (and selfintersections) situation becomes 



more complicated [TT[]. Namely, each cusp leads to its own local Z - fac- 
tor, depending on the cusp angle (but not on another characteristics of the 
contour). The Eq.(p!5D for the loop with K cusps having angles 7^ is to be 
modified in the following way: 

K 

W(C) =n^.(70-W^ren(C^) (17) 
i=l 

The mixing between linear and logarithmic divergencies, mentioned above 
has some specific features for the loops with cusps. Namely, each cusp in ad- 
dition to the linearly divergent factors, removed by the test particle mass 
renormalization introduces logarithmically divergent terms which have noth- 
ing to do with the coupling constant renormalization and must be subtracted 
by their own Z(7j)-factors. There are two types of terms of this kind, which 



do depend on the cusp angle and which do not [Til]. '^^^ latter plays im- 
portant role in the perturbative expansion of the correlators. To see this, 
note, that the correlators can be easily obtained from the Wilson loop func- 
tional. Indeed, differentiating (|I|) with respect to bOy,y[xi) and assuming that 
the contour C connects the points m, v with the reference point Xq along the 
straight lines, one finds for the bilocal correlator: 

^Tr/((F,.(..,Xo)F,.(t;,Xo))) = -j/^^^j-. (18) 

and analogous expressions for higher cumulants. Note, that area of the min- 
imal surface bound by the contour C from (0) is equal to zero and this 
contour generally has at least two cusps (typically even four, if the reference 
point Xq does not coincide with one of the correlator arguments). 

Comparing Eqs. (|l^) and (p!8|) one concludes, that the derivatives in the 
r.h.s. of (|l^) must act on renormalized part of the Wilson loop as well as 
on the Z - factors. As a result mixing between different terms appears and 
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contour divergencies cannot be separated in some simple way as they can in 
the case of loops, so the perturbative behaviour of the correlators is controlled 
by both ultraviolet and contour logarithmic divergencies. 

Let us do the final remark. It may be easily shown, that taking into 
account only the bilocal correlator in the cumulant expansion (|^) one has 
perturbative contributions to the string tension at the one loop level since 
the kroneker part of (||) does not equal to zero. But area law is perturba- 
tively impossible. The problem is cured by the contributions from the higher 
cumulants, in particular triple correlator (FFF) also contain such kroneker 
parts being considered at the same {g"^ in our case) order of perturbation 



theory |T^. As a result string tension is equal to zero in the field-strength 
formulation as it must be at any given order of perturbation theory if all cu- 
mulants contributing to string tension at this order are taken into account. 
Note in this respect that there are no terms proportional to the /3-function 
coefficients in the function D 0] - it is a direct consequence of the fact that 
nonzero D has no physical sense in perturbation theory and all contributions 
into it must be cancelled by the next-order terms. 
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